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Abstract. We have obtained a sample of 143 normal galaxies with X-ray luminosity in the range 10 38 
s _1 from the cross-correlation of the ROSAT HRI Brera Multi-scale Wavelet (BMW - HRI) Catalogue with the 
Lyon-Meudon Extragalactic Database (LED A). We find that the average X-ray properties of this sample are in 
good agreement with those of other samples of galaxies in the literature. We have selected a complete flux limited 
serendipitous sample of 32 galaxies from which we have derived the logN-logS distribution of normal galaxies 
in the flux range 1.1 — 110 x 10~ 14 erg cm -2 s _1 . The resulting distribution is consistent with the euclidean 
— 1.5 slope. Comparisons with other samples, such as the Extended Medium Sensitivity Survey, the ROSAT All 
Sky Survey, the XMM - Newton/2dF survey and the Chandra Deep Field Survey indicate that the logN-logS 
distribution of normal galaxies is consistent with an euclidean slope over a flux range of about 6 decades. 

Key words. X-ray - galaxies: general - X-ray surveys 



1. Introduction 

Detailed X-ray studies of normal galaxies have been 
possible only with the advent of imaging instruments 
aboard the Einstein observatory. The Einstein results 
we re summarized i n the catalogue and atlas published 
bv IFabbiano et all 1^992), who constructed a large ho- 
mogeneous sample of 493 galaxies, included in cither 
A Revised Shapley-Ames C atalog of Bright Galaxies 
ijSandaee fc TammannL Il98llb or in the Second Revise d 
Catalog of Bright Galaxies l)de Vaucouleurs et al.l fl976) . 
both targets and serendipitously detected in Einstein 
field s. While representative of the galaxy pop ulation (see 
e.g. lShaplev et alll200ltlEskridge etailll995[) . it was not 
constructed to be a complete, unbiased sample, and it is 
likely to contain unknown selection biases. 

Elliptical galaxies were found to retain large amounts 
(10 s - 10 11 M ) of hot gas (T - 10 7 K) whose thermal 
emission dominates their X-ray luminosities, while in nor- 
mal spirals the integrated contribution of the evolved stel- 
lar sources, such as supernova remnants and X-ray bina- 
ries, i s generally the dominant component (se e iFabbianol 
Il989t IFabbiano et"afl Il992i iKim et all Il992fl . Extended 
emission from a hot, gaseous com ponent in spiral galaxies 
was detected only in some cases l|Fabbiano fe Trinchieri 



Il987t IVogler fc Pietschl Il996t ITrinchieri et"all Il988|h or 

associated with starburst activity. 

Subsequent observations of individual sources made 
by ROSAT and ASCA confirmed most of the Einstein 
results and added interesting information on the 
X-ray properties of normal ga l axies (see among 
others iRoberts fc WarwickL l2000t iRead et all Il997t 
iBrown fc Bregmanl Il998t iBeuing et all Il999|) "in the lo- 
cal universe. 

With the launch of XMM - Newton and Chandra, 
the study of the X-ray properties of "normal" galax- 
ies at interme d iate ( z ~ 0.1) or cosmolog i cal di st ances 
llBrandt et all 1200 it iHornschemeier et all 12002 . l2003t 
iGeorgakakis et all l2003t l2004albt INorman et ail l2004 



was made possible, thanks to significantly improved sen- 
sitivity, spatial and spectral resolution of the instruments. 
In spite of the large number of papers, however, a truly 
complete sample of X-ray emitting normal galaxies in 
the local universe with a significant number of objects 
has not been properly d iscussed in the literature so far. 

galaxy number density at 
1 from stacking analysis, 
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F n 5- Q ~ 7 X 10- 

while lGeorgakakis et al.l l|2004af) present a sample of 26 X- 
ray sources detected in an area of ~ 2.5 deg -2 , of which 
only 2 are however classified as n ormal galaxies. Only re- 
cently IGeorgakakis et all l)2004bj) have presented a larger 
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sample of 11 normal galaxies detected in an area of ~ 4.5 
deg~ 2 . 

Two other samples of "normal" galaxies are avail- 
able in the literature, selected in the Chandra 
Deep Fields l|Hornschemeier et all 120031: iNorman et all 
2004). However the i r med ian redshifts (z = 0.297, 
Hornschcmeicr et al. and z = to z — 1.3, 

Norman et all (|2004j)) indicates that they should not be 
considered as "local". 

The large database provided by ROSAT has been ex- 
ploited only margina lly to derive unbias e d and complete 
sample of galaxies. IZimmermann et "ahl ()200l[) have se- 
lected a sample of candidate normal galaxies from the 
ROSAT All Sky Su rvey (RASS) Bright Source Catalogue 
i Voges et al!ll999l) above a flux limit of about 10~ 12 erg 
cm~ 2 s _1 (0.1 — 2.4 keV band). A few samples have been 
derived from the ROSAT Position Sensitive Proportional 
Counter (PSPC) pointed observatio ns such as the WGA 
(White. C T iommi. Rr. AngelTnl Il994h and the ROSPSPC 
(ROSAT team 2001), and a new catalogue of galaxies is 
in progress (G. Peres, private communication). Here we 
exploit the potential provided by the Brera Multi-scale 
Wavelet (BMW - HRI) catalogue l|Panzera et all [2003) 
to extract a sample of normal galaxies, as we discuss in the 
next sections. While the PSPC is probably more efficient 
at detecting faint and extended sources such as galaxies, 
the sharp core of the HRI point spread function allows the 
detection of sources in more crowded fields and to estab- 
lish the extension for bright small-size sources, providing 
a good complement to the PSPC data. 



2. The sample 

In order to create a complete, serendipitous sample of 
galaxies with X-ray emission, we made use of X-ray 
data from the BMW ROSAT HRI catalogue and optical 
data from LEDA (Lyon-Meudon Extragalactic Database). 
The BMW - HRI catalogue consists of 29089 X-ray 
sources detected in 4303 ROSAT HRI pointed fields with 
exposure times lo nger than 100 s us i ng a multiscale 
wave l et algorithm (Lazzati et all Il999t ICamoana et all 
Il999l IPanzera et alll2003|) . Sources detected with a sig- 
nificance > 4.2 a are contained in the catalogue, that pro- 
vides name, position, count rate, flux and extension along 
with the relative errors. In our study we used the full cat- 
alogue, but we excluded X-ray sources in the Trapezium 
field, which is a rich stellar cluster in the Milky Way, where 
the high density of the X-ray sources would prevent proper 
optical identifications. The BMW - HRI catalogue can 
be searched via the HEASARC Browse 1 or via the Brera 
Observatory web site. 2 



1 http : //heasarc . gsfc.nasa.gov/cgi-bin/W3Browse/ 
w3browse .pi 

2 http : //www.merate .mi . astro . it : 8081/interroga/ 
dbServer?cmd=bmw2 



Created in 1983 at Lyon Observatory, LEDA 3 has been 
the first database of extragalactic objects and it is contin- 
uously updated. It gives a free access to the main astro- 
physical parameters (coordinates, morphological type, di- 
ameter and axis ratio, apparent magnitudes and colors, ra- 
dial velocity, surface brightness, etc) for about 10 6 galax- 
ies over the whole sky. The completeness i n apparent B - 
magn itude is satisfied up to = 15.5 fsee lPaturel et all 
ll 997ft . 

To obtain a representative sample of galaxies we 
started from the BMW - HRI catalogue and included only 
serendipitous detections, avoiding the targets. We chose a 
3' radius to define the typical region of the target, and 
selected only sources at off-axis angles 9 > 3'. In spite 
of their off-axis location, we had to exclude 12 additional 
sources that were targets of the observations. We then 
cross - correlated the positions of the X-ray sources in 
the BMW - HRI catalogue with those of galaxies present 
in the LEDA database version of 1999, with a tolerance 
of 20", which should be a reasonable guess to detect ex- 
tended objects like galaxies and to avoid most of chance 
coincidences. This criterium is not appropriate for very 
extended galaxies (like e.g. M 31 or M 33), where a large 
number of sources is detected at distances significantly 
larger than our tolerance radius. Therefore we could be 
selecting against large size galaxies, if there is no source 
within a distance of ~ 20" from the nucleus. However, 
this should not be a concern in this study, since the sur- 
face density of large galaxies is small: in particular, in the 
LEDA catalog, the density of galaxies with D 2 5 > 3' is 
< 3 x 10 -2 deg -2 , which implies about < 10 in the area 
we surveyed. 

To check the goodness of our choice, we plot in Fig. 
^the relative shifts between optical and X-ray positions. 
We use only point-like X-ray sources because in extended 
sources the association with a single optical object could 
be misleading (e.g. in groups and pairs). About 90% of 
identifications are within 13" (represented by the dashed 
circle in Fig. QJ. 

This is in good agreement with the HRI po- 
sitional uncertainty, since the best attitude solu- 
tion guarantees that on average known objects are 
detected within 10" of their catalog position, al- 
though with possible additional discrepancy mainly 
in declination (see the ROSAT Handbook e.g. at 
http : //heasarc . gsf c .nasa.gov/docs/rosat/ruh/ 
handbook/handbook. html). We verified that all X-ray 
sources with X-ray/optical off-set > 10" are still within 
the galaxy (i.e. D25). We conclude that the associa- 
tion of the X-ray sources with galaxies found by the 
cross-correlation is sound on positional grounds. 

The cross-correlation yields 399 X-ray sources asso- 
ciated with 281 galaxies. Being constructed field- by-field, 
the BMW - HRI catalogue contains multiple detections of 
the same source. Therefore we only have 283 distinct X- 

3 There exist several mirrors to access LEDA; we used the 
OAB one, http://www.brera.mi.astro.it/hypercat/ 
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Fig. 1. Displacement between optical and X-ray positions 
for point-like sources in the total sample. The solid circle 
refers to the cross - correlation radius of 20"; the dashed 
one, with 13" radius, contains about 90% of the identifi- 
cations. 



ray sources associated with galaxies (two distinct sources 
but of different extent are associated with both NGC 1399 
and M 86; we list them in Table but we consider only 
the largest one for computing fluxes and luminosities). 

To check the results of the cross-correlation, we in- 
spected HRI images for each source and obtained informa- 
tion from the NASA/IPAC Extragalactic Database (NED) 
to a) eliminate spurious coincidences [15 sources], b) elim- 
inate AGNs [47 objects] and c) select clusters [57 objects]. 

a . We eliminated a source when the obvious optical 
counterpart was not a LEDA galaxy, but a background 
or foreground object. This was checked both using 
NED and the X-ray/optical contour maps. 

b We excluded object s classified in 

IVeron-Cettv fc Veron Catalogue! l)200l|) as Seyfert 1 
galaxies, QSOs, BL Lacs or AGNs. For Seyfert 2 
galaxies, the X-ray emission could have a non nuclear 
origin, so we eliminated them when we saw from 
images that the X-ray emission was point-like and 
well centered on the nucleus, while we retained those 
for which we found an extended emission. 
We also retained starburst galaxies, LINERs and 
objects that we know from the literature could have a 
nuclear source but also non - nuclear X-ray emission 
(e.g. NGC 3079). 

c . Several galaxies lie inside a cluster: we inspected 
X-ray images and we eliminated those objects whose 
X-ray emission was indistinguishable from the clus- 
ter's, but we retained galaxies for which emission 



clearly associated with the galaxy is detected above 
the cluster background (see e.g. PGC 12350 in Fig.|2J). 
We retained galaxies in poor groups even when 
they are the brightest member since there is 
still ambiguity in the literature between emission 
from bright early- type galaxies and groups, which 
are often analyzed in the same context (see e.g 
Mamorl Il992t IPelF Antonio et all 11994 IPildis et al 



199.51 IPonman et all Il996t iMahdavi et all Il997: 
Mulchaev fc Zabludofll 1^; iHelsdon fc Pontnan. 
2000 iMulchaev et all 2003t lOsmond fc Ponman . 
2004UHelsdon fc Ponmanl l2003t Ijones et alll2003|l 



Moreover, in 5 cases the same X-ray source was asso- 
ciated with two or more galaxies in a pair or in a group: 
since we could not discriminate on a positional basis, we 
chose the brightest galaxy in the pair or in the group. 

These selection criteria yielded a total of 195 X-ray 
sources (including multiple detections) associated with 
143 galaxies whose properties are listed in Tabled avail- 
able in electronic form. We report here the first page as 
an example. Column (1) gives the BMW - HRI name of 
the source, columns (2) and (3) the position of the X- 
ray peak, column (4) labels the source as point-like (p) 
or extended (e), column (5) gives the radius in which we 
computed count rates for extended sources (see section|3J), 
column (6) the major semiaxis of the ellipse, for extended 
X-ray sources for which count rates have been computed 
in elliptical regions (see section yj, column (7) the minor 
semiaxis of the ellipse, column (8) the LEDA galaxy as- 
sociated to the X-ray source, column (9) the name of the 
galaxy in other common catalogues (e.g. NGC), column 
(10) the morphological type, column (11) the distance of 
galaxy, column (12) the apparent B - magnitude corrected 
for galactic extincti on, inclination and redshift effects (see 
IPaturel et ailll997|) . columns (13) and (14) the HRI count 
rate and error, columns (15) and (16) the X-ray flux and 
error (0.1 — 2 keV; count rates and fluxes have been re- 
computed with respect to values reported in the BMW - 
HRI catalogue; see section for details), column (17) the 
logarithm of X-ray luminosity and column (18) specifies 
whether the galaxy is in the complete subsample (c), in 
cluster (CI) or in group (Gr). 

Information about magnitudes and redshifts are from 
LEDA and NED except for 3 objects that we observed 
ourselves (see section EJ. We calculated distances from 
redshifts, assuming Hq — 50 km s _1 Mpc -1 ; when the 
heliocentric radial velocity of galaxy was less than 3000 
km s , ho wever, we us ed distances from Nearby Galaxy 
Catalogue (jTullvL Il988h corrected for H Q = 50 km s" 1 
Mpc" 1 . 

One hundred and sixteen of the 143 galaxies have 
redshift and magnitude (including our observations). 
Nineteen galaxies have mg and no z, 3 have z and no 
mj and 5 have neither z nor m b ■ The reshift distribution 
ranges from z = from z ~ 0.15, but ~ 90% of galaxies 
have z < 0.07. 
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(1) BMW - HRI source name H 

(2) X-ray coordinates: RA (J2000) |- 

(3) X-ray coordinates: Dec (J2000) o 

(4) X-ray source point-like (p) or extended (e) p_ 

(5) Radius of circle in which count rates were computed, for extended X-ray sources ^ 

(6) Major scmiaxis of ellipse in which count rates were computed, for some peculiar extended X-ray sources » 

(7) Minor semiaxis of ellipse in which count rates were computed, for some peculiar extended X-ray sources 

(8) LEDA galaxy associated ™ 

(9) Other galaxy name ^ 

(10) Morphological galaxy type g 

(11) Galaxy distance (Mpc) ^ 

(12) Apparent B - magnitude 

(13) X-ray count rate (10~ 3 count s _1 ) pi 

(14) X-ray count rate error (10~ 3 count s _1 ) O 

(15) X-ray flux (0.1 - 2 kcV band; 10~ 13 erg cm" 2 s" 1 ) | 

(16) X-ray flux error (0.1 - 2 kcV band; 10~ 13 erg cm" 2 s -1 ) |f 

(17) Logarithm of X-ray luminosity (0.1 — 2 keV band; erg s _1 ) a 

(18) Notes: c if the galaxy is in the complete subsample, CI if in cluster, Gr if in group 



Table 1. First page of Table of the total sample, available in electronic form. 
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2.1. The complete serendipitous subsample 

To study the general properties of the sample we need 
to derive a subsample with well known completeness cri- 
teria and limits. To this end we constructed a complete 
sample with both X-ray and optical flux limits. We must 
consider both the X-ray and the optical completeness cri- 
teria. The X-ray completeness is related to the BMW - 
HRI catalogue, that includes all sources with a significance 
> 4.2 a. To take into account the optical limits, we ex- 
cluded galaxies fainter than m B = 15.5 (th e completeness 
limit assumed for the LEDA Catalogue, see lPaturel et all 
Il997h . 

We have also excluded objects at low galactic latitude 
< 10°) to avoid source confusion in the galactic plane 
and X-ray sources with off- axis angle 9 > 18' to match 
the circular HRI field of view used in the sky coverage 
computation (see section 0). 

The resulting sample of 96 objects is complete both 
in X-ray and optical, at the given limits. However, since 
galaxies are often in agglomerates, some of them are ex- 
pected to be related to the targets and therefore not truly 
serendipitous. We therefore excluded all sources known to 
be associated with the target (e.g. galaxies in pairs, groups 
or clusters; 52 objects). When an association was not doc- 
umented (e.g. from NED, LEDA) we conservatively ex- 
cluded galaxies at the same redshift of the target (12 ob- 
jects). 

The complete, serendipitous sample of 32 galaxies thus 
obtained is given in Table |3| and will be used to calculate 
the logN-logS distribution in the local universe (z < 0.07), 
as will be described in detail in section 



3. X-ray characterization 

The BMW - HRI catalogue provides count rates derived 
using the wavelet algorithm, in an automat ed way and un- 
der particular assumptions (for details, see lLazzati et all 
ll999tlCampana et allll999l) . We verified that, while count 
rates are correctly computed for point- like sources, for 
very extended sources such as NGC 1399, the extension 
and count rates given by the algorithm are underesti- 
mated. Moreover, when multiple observations are available 
we can improve the statistics by considering the full set 
of data. We have therefore recalculated all count rates, 
using the BMW - HRI positions and original HRI data 
retrieved from the public archives 4 , in pulse height an- 
alyzer (PHA) channels 1 — 10 to increase the signal-to- 
noise ratio (for a just ification of this choice of PHA, see 
iTrinchieri et allll997|) . When multiple observations of the 
same field are available, we summed the data if they have 
the same pointing coordinates and comparable exposure 
times. Otherwise, we typically used the longer exposures 
or those where the source is closer to the field's center. 



http://wave.xray.mpe.mpg.de/ and 
: //heasarc . gsf c . nasa. gov/ 



3.1. Count rates 

We classified sources as point-like or extended, based on 
the radial distribution of the emission relative to the shape 
of the HRI point sprea d function (P SF; for a description 
of the ROSAT PSF see lBoesell2000|) at the corresponding 
off-axis angle. For the "extended sources" the counts are 
taken from the largest region that contains source counts, 
generally a circle of radius reported in Tabled evaluated 
from the radial profile. In few cases, given the particular 
shape of the source, we used ellipses to evaluate the counts; 
we give in Tabled minor and major semiaxis. There are 
also a few extended sources for which ad hoc regions (nei- 
ther circles nor ellipses) have been used to estimate counts: 
these are PGC 5323, PGC 5324, PGC 13418, PGC 13433 
and PGC 56962. 

For point-like sources the counts are obtained in a cir- 
cular region centered at the peak of X-ray emission, with 
radius that includes about 90% of source counts accord- 
ing to the PSF. The PSF degrades as the angular distance 
from the center of field increases, so we chose a radius of 
r = 18" for 3'< 6 < 10', r = 25" for 10 ' < 9 < 15' and r = 
40" for 9 > 15', following iBoesel (|2000j) . We evaluated the 
background in an annulus concentric to the source radius, 
with radii depending on the off-axis. When the source was 
particularly faint, we calculated count rates in a circle of 
radius corresponding to a smaller fraction of the PSF, to 
increase the signal-to-noise ratio . We then co rrected the 
count rate accordingly, following Bocsc (2000). 

We compared count rates obtained in this way with 
those reported in the BMW - HRI catalogue and we found 
a general agreement, with the exception of some sources 
whose extension had been largely underestimated by the 
wavelet algorithm, as stated above. Indeed, our count rates 
would be the equivalent to the "counted count rates" re- 
ported in the catalogue, rather than those computed with 
the wavelet algorithm. In the comparison with this quan- 
tity we found a systematically higher count rate, consis- 
tent with the larger PHA interval used (1 — 10 in our 
analysis and 2 — 9 in the BMW - HRI catalogue). 

The resulting net count rates are given in Table 
corrected for vignetting and lost counts due to the PSF 
(for point sources only). 

3.2. Fluxes and luminosities 

The count rates were converted into 0.1 — 2 keV fluxes us- 
ing a conversion factor corresponding to a bremsstrahlung 
spectrum with kT = 5 keV plus the line of sight absorp- 
tion a ppropriate for each source from lDickev fc Lockmanl 
l|l990l) (reported in Table EJ. Although this spectrum 
might not be suitable for all kinds of sources, the flux 
in the ROSAT energy window depends only weakly on 
the spectral model assumed, while it is more dependent 
on low energy absorption. The resulting fluxes are in the 
range 10~ 14 — 10~ n erg cm~ 2 s _1 . 



http 
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N R (cm-) 


CF (erg cm 2 count 1 ) 


1 x 10 2U 


3.7 x 1(T U 


2 x 10 20 


4.1 x KT 11 


5 x 10 20 


4.9 x KT 11 


8 x 10 20 


5.5 x KT 11 


3 x 10 21 


9.3 x 10" 11 



Table 2. Galactic N-a and corresponding unabsorbed flux 
in the band 0.1-2 keV for 1 count/s, assuming a thermal 
bremsstrahlung spectral model with kT = 5 keV. 



To calculate luminosities, we used distances listed in 
Tabled The corresponding range in Lx is 10 38 — 10 43 erg 



4. The atlas 

We provide in Fig. an overlay of X-ray contours of the 
detected galaxies onto optical images from the Digital Sky 
Survey II (DSS II) available from the ESO archive 5 . The 
Figure is available in electronic edition. We report here the 
first page as an example. When available, we used optical 
images obtained with the blue filter, otherwise we used 
those obtained in the red filter. For PGC 209730 only the 
DSS I plate is available. 

Galaxies, ordered in RA and generally at the center of 
the field, are identified by their PGC name. 

X-ray contours are produced from images in the PHA 
range 1 — 10, smoothed with a Gaussian function with a 
= 5" for point-like sources and with a = 10" for extended 
sources. 



5. Optical observations 

In order to measure redshifts and magnitudes for some of 
the galaxies of our sample, we made spectroscopical and 
photometrical observations at the 1.52 meter telescope 
of the Osservatorio Astronomico di Bologna, at Loiano 
(Italy) on the nights of the 16th and 17th October 2001. 
Because of bad atmosferic conditions, we were able to ob- 
serve only 3 galaxies. A spectrophotomctric calibration 
star was also observed. We present the results obtained in 
Appendix lAl 



6. Comparison with literature results 

To verify whether the total sample of 143 objects is rep- 
resentative of the X-ray properties of normal galaxies, 
we calculated X-ray luminosities where possible (Table 
^) and we plotted the distributions of Lx , Lb, the ratio 
Lx/Ls and the Lx-Lb relationship for spiral and early- 
type galaxies. For the 19 galaxies for which the redshift is 
not known but we have mg, the ratio is calculated from 



40 - 



30 - 



20 



10 




log(F x /F B ) 

Fig. 3. X-ray - to - optical ratio distribution for the total 
sample. Both fluxes are in units of erg cm~ 2 s . 



fluxes. We compared our results with those in the litera- 
ture and found a good general agreement. 
In particular: 



1. the bulk of the galaxies in the sample has an X-ray lu- 
minosity between 10 3 8 erg s 
accordance with, e.g.. lFabbianol l|l98 



1 and few 10 42 erg s 1 



m 

7 objects have 
Lx ^ 10 43 erg s _1 , but most of them lie in a group, 
so that the intcrgalactic medium could contribute to 
their luminosity; we cannot exclude the presence of an 
unidentified AGN for some of these objects. 
For spiral galaxies we found a linear relationship be- 
tween X-ray and o ptical luminositie s (Ly oc L^ 0±0 ' 2 ), 
in agreement with iFabbiano et alJ l|l992() . In a sub- 
sequent, more complete statistical analysis of the 
23 4 "normal" spir a l and i rregular galax i es rep orted 
bv IFabbiano et alJ |l992l). IShaplev et all (|200l|) and 
IFabbiano fc Shaplevl (2002) find however that the Lx 
— Lb relationship is significantly steeper than linear, 
with a slope of about 1.5. Our study is based on a much 
smaller sample (32 spiral and irregular galaxies), and 
small statistic s could account for the disc repancy, as 
pointed out bv IFabbiano & S haolc 1 ll2002h m compar- 
ing their findings with previous works. 
For early-type galaxies we found a steeper relation- 
shi p, Lyoc Lb 1 - 6 ^ - 2 , co nsist ent with those o btain ed 
bv IFabbiano et alJ (|l992j) and lEskridge et alJ l)l995|) . 
Values of the X-ray - to optical ratio in our sample 
cover roughly the same range (—4 < \og(Lx /Lb ) < 
with Lx and Lr in erg s" 1 ) as the spiral galaxies 



http : // archive . eso . org/dss/dss 



inlShaplev et all (l2001fl and the early-type galaxies in 
iFskrinVe et al] lll995h . but their distributions are dif- 
ferent. We plot in Fig.Ejthe histogram of the X-ray - to 
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Fig. 4. Sky coverage for the fields at galactic latitude |6| > 
10° and off-axis angle 3'< 9 < 18', computed by assuming 
a thermal bremsstrahlung spectrum with kT = 5 keV and 
galactic line of sight absorption. 




lo g S o.5-2.o ( er § cm " S s_1 ) 



Fig. 5. The integral logN-logS distribution for the com- 
plete serendipitous subsample (asterisks); the solid line 
represents the —1.1 slope and the dashed line the euclidean 
— 1.5 slope. 



- optical ratio of the total sample of 143 galaxies. We 
will analyze this subject in greater detail in subsection 

rm 

7. The logN-logS distribution 

The complete serendipitous sample of galaxies derived in 
section ( listed in Tabled was used to calculate the in- 
tegral flux distribution (logN-logS) of normal galaxies with 
X-ray emission above the X-ray flux limit of the BMW - 
HRI Catalogue and B - magnitude < 15.5. 

The sensitivity of the HRI instrument is not uniform 
over the entire field of view. Moreover the observing time 
is different for different fields so we need to calculate the 
area surveyed at any given flux (sky coverage) . 

In the BMW - HRI catalogue the published sky 
coverage was calcul ated by means of simulations (see 
IPanzera et all l2003|) . In this work we used a sky cover- 
age calculated with the same procedure, but with param- 
eters that reflect our source selection criteria. Therefore 
we included only fields with galactic latitude \b\ > 10°, 
we considered only an annular region with 3'< 9 < 18' 
(the lower limit is to account for the target region; the 
upper limit is the largest radius within the field of view of 
the detector in the assumption of circular symmetry) and 
we assumed a bremsstrahlung spectrum with kT = 5 keV 
plus the line of sight absorption. 

The resulting sky coverage is plotted in Fig. H The 
maximum area is ~ 314 deg 2 and corresponds to fluxes 
above ~ 10 -12 erg cm' 2 s _1 . The surveyed area is ~ 196 
deg 2 at 10~ 13 erg cm~ 2 s _1 and ~ 3 deg 2 at ~ 1.1 x 



10 -14 erg cm' 2 s' 1 (the lowest flux for the galaxies in 
our sample). 

For consistency with the sky coverage calculation, the 
fluxes given in Table |31 are derived from the original count 
rates estimated in the BMW - HRI catalogue using the 
wavelet algorithm and PHA channels of HRI from 2 to 9, 
which are corre c ted fo r vignetting and PSF according to 
ICamnana et al.l 1 1999) and computed for the energy band 
0.5 - 2 keV, the standard of BMW - HRI Catalogue. 

The integral logN-logS distribution of the sample is 
shown in Fig.[S]and covers two decades in flux, from ~ 1.1 
to ~ 110 x 10~ 14 erg cm~ 2 s _1 . The overall distribution 
could be approximated with a slope of ~ —1.1 (solid line 
in Fig. EJ). However the euclidean slope of —1.5 (dashed 
line) is also consistent with the data: the excess of galaxies 
at the highest fluxes is small and consistent within the 
limited statistics. Moreover, at the lower fluxes, we could 
have some problems with incompleteness, as we discuss 
below. 



7.1. Comparison with the literature: X-ray 
optical ratio distribution 



to 



Before comparing our logN-logS distribution with those 
of other samples in the literature, we need to better in- 
vestigate the X-ray - to - optical ratio distribution of 
our complete sample. We will consider samples derived 
from ROSAT and Einstein observations that cover a flux 
range similar to ours. We will include samples derived from 
XMM - Newton and Chandra surveys, that cover a flux 
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LEDA name 
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Table 3. Galactic column density (from Dickey & Lockman 1990), BMW count rates and fluxes for galaxies in the 
complete serendipitous sample, ordered by increasing flux. 



range significantly fainter than ours, in the discussion of 
the logN-logS (subsection E3- 

The best available comparison could be with 
the sample of candidat e normal galaxies found by 
IZimmermann et all l|200lj) in the ROSAT All Sky Survey 
(RASS) and with the normal galaxies found in the 
Einstein Extended Medium Sensitivity Survey (EMSS; 
iGioia et "all Il99d(l . both X-ray selected. An effectively 



galaxy sample i 


7 abbiano et alJ 


■ Il992: Shaolev et al.Ll2001: 


Eskridee et all 


1995). 



IZimmermann et alJ ()200l[) made a correlation study 
of the RASS Bright Source Catalogue llVoees et allll99al 
with the Catalogue of Principal Galaxies IjPaturel et all 
from which they selected a sample of 198 candidate 
galaxies, i.e. X-ray sources whose optical counterpart was 
not designated as AGN in the literature. These selection 
criteria are similar to ours, and the Catalogue of Principal 
Galaxies is a preliminary version of the current LEDA 
database, so the two samples can be easily compared; how- 



ever most of the IZimmermann et"al1 ((2001) sources have 
fluxes above 10~ 12 erg cm~ 2 s _1 (computed in the 0.1 — 2.4 
keV band, assuming a power law spectrum with photon 
index T = 2.3). 

The EMSS was obtained from the analysis of 1453 im- 
ages of the imaging proportional counter (IPC) on board 
the Einstein Observatory. The survey covers an area of 
778 deg 2 at \b\ > 20° with limiting sensitivity ranging 



from ~ 5 x 10 



to 



(0.3-3.5 



keV band). 835 serendipi tous sources were detected at or 
abov e the 4 a level (see IGioia et "all Il99f)t IStocke et all 
Il99lh . Among these, 17 were identified as normal galaxies. 

The Einstein sample is the catalogue of normal 
galaxies observ e d by the Einstein satellit e, compiled by 
iFabbiano et alJ l|l992|) and reanalyzed by lE^kridge^et^l' 
(1995, for early-type galaxies) and bv IShaplev et al 
(200lj, for spiral galaxies); we no longer distinguish here 
the two samples since early and late type galaxies are 
mixed in the BMW - HRI sample. 
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Fig. 6. B flux versus X-ray flux for our total sample (solid 
triangles), for the Einstein sample (empty circles), for the 
Zimmermann et al. 2001 sample galaxies (asterisks), and 
for the EMSS galaxies (empty pentagons). Solid lines cor- 
respond to log(Fx/Fs)= 0, —2, —4, as indicated; the hor- 
izontal dashed line is the nig = 15.5 limit. We have consid- 
ered only detections in the Einstein sample: the exclusion 
of upper limits should not affect the results, since their 
distribution is consistent with that of the detections (see 
Fig. 5 in Shapley et al. 2001). We do not plot objects at 
Fo.1-2 < 2 x 10~ 12 erg cm~ 2 s _1 in the Zimmermann et 
al. (2001) sample because below this flux their sample is 
not complete (see text). 

In Fig. we plot the distribution of Fx and Fb val- 
ues from all samples considered. For consistency with the 
values in Tabled all fluxes are converted to the 0.1 — 2 
keV energy band, using a thermal bremsstrahlung spec- 
trum with temperature kT = 5 keV and N# = 3 x 10 20 
cm . 

Fig. [f)| clearly indicates that galaxies belonging to dif- 
ferent samples populate different regions in the plot. Our 
sample (solid triangles) has an X-ray flux range between 
10 -14 and 10 -11 erg cm~ 2 s _1 and log(Fx/Fs) between 
—4 and (see also Fig- EJt - The bulk of Einstein galaxies 
(empty circles) is tipically at higher average fluxes and at 
lower values of log(Fx/Fs) (between —4 and —2) com- 
pared to our distribution and the X-ray selected samples 
in general. This sample is the largest, it is effectively op- 
tically selec ted and reasonably c lean of the contamination 
from AGN l|Shanlev et alll200lh . However, since it is not 
complete, it might not provide the true distribution of the 
X-ray - to - optical ratios. 

The second la rgest sample is derived from 
IZimmermann et"al1 l|200lh . with an additional flux 
limit Fx > 2 x 10~ 12 erg cm~ 2 s _1 (for completeness, 



see discussion below) and an upper limit in luminosity 
at logLx = 42.7 erg s _1 (the highest luminosity in our 
complete sample), comparable to their limit to exclude 
potential AGNs fr om the sample. The d i stribu tion of can- 
didate galaxies of IZimmermann et*aT] l|200lh (asterisks) 
is significantly different from that of the Einstein sample 
and extends at log(Fx/Fs) > 0. 

The Fx/Fb distribution for EMSS galaxies (empty 
pentagons) is at intermediate values and more consistent 
with that of our sample. 

Also plotted in Figgis the optical flux limit applied 
to our sample, nig = 15.5. It is evident that the exclusion 
of galaxies fainter than this for the logN-logS calculation 
has an effect that increases as flux decreases. We have at- 
tepted to quantify it in order to correct the curve for lost 
objects. Unfortunately we could not properly estimate the 
correction because none of the galaxy samples available al- 
lows us to derive the true distribution of X-ray - to - opti- 
cal ratios of normal galaxies. The EMSS should represent 
the true Fx/Fb distribution, but, given the small num- 
ber of galaxies (17, of which only 15 with B-magnitude), 
statistical errors are large. The difference in the distribu- 
tion of ratios between the two larger samples in Fig. H3 
sug gests that they might be affected by opposite biases: 
the IZimmermann et all l|200lh sample is likely to contain 
unclassified AGNs, while the Einstein sample could lose 
objects at the highest X-ray - to - optical ratios. If we use 
the three samples to estimate how many galaxies are lost 
as a function of X-ray flux, we find that the corrections 
to the logN-logS are small and the recomputed curve is 
consistent with the euclidean slope in the observed flux 
range. 

We also tried to estimate the correction by considering 
optically fainter galaxies. At mp = 16, LED A is about 
90% complete (see Fig. 7 in IPaturel et all Il997h : if we 
include galaxies down to this flux limit, we only add 4 
objects to our serendipitous sample, distributed over the 
whole range of X-ray fluxes. Their inclusion only has a 
marginal influence on the normalization, and none on the 
slope of the distribution. 

We conclude that, since the effects introduced by the 
optical limit are small, the logN-logS we derive is consis- 
tent with the euclidean slope. 

7.2. Comparison with the literature: logN-logS 

Figure [3 shows the comparison between the logN-logS de- 
rived above with several available from the literature. All 
fluxes are recomputed in the 0.5 — 2 keV range that we 
use. 

We find an excellent agreement with other samples 
that cover similar or brighter flux ranges than the present 
sample. 

The candidate galaxies ijZimmermann et alll200lf) ap- 
pear to connect smoothly with the euclidean extrapo- 
lation of the BMW - HRI logN-logS above 10~ 12 erg 
cm~ 2 s _1 . We interpret the flattening observed in the 
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Fig. 7. Comparison of the logN-logs distribution obtained 
from the present sample (asterisks) with several other re- 
lations for normal galaxies from the literature. The dashed 
line respresents the euclidean relationship of Figure [3] ex- 
trapolated to the full flux range considered. Symbols re- 
fer to: EMSS galaxies (empty circles); candidate galax- 
ies in Zimmermann et al. 2001 (solid squares); nor- 
mal galaxies (crosses), including LLAGN (solid circles) 
(Georgakakis et al. 2004a), and stacking analysis (solid 
triangles, Georgakakis et al. 2003), found in the XMM- 
Newton /2dF survey; XMM-Newton /Sloan DSS normal 
galaxies from Georgakakis et al. (2004b) (empty squares); 
Chandra Deep Fields results from Hornschemeier et al. 
(2003) (empty triangles); optimistic or pessimistic curves 
given by Bauer et al. (2004) (upper/lower dotted line); 
Bayesian sample of Norman et al. (2004) reported by 
Ranalli et al. (2004) (dot - dashed line). Open stars and 
horned-shaped box indicate the detections and fluctuation 
analysis predictions of Miyaji & Griffiths (2002). 



IZimmermann etlrtl l|200l[) data at lower fluxes as a re- 
sult of their selection criteria. In fact, they indicate a > 
90% completeness for count rates > 0.1 count s _1 , which 
converts to a flux of Fo.5-2 ~7x 10~ 13 erg cm~ 2 s -1 . 

Although small, the EMSS is a truly complete sam- 
ple, since it is sercndipito usly X-ray selected and virtu- 
ally completely identifi ed ( Gioi a et al l ll99f]UStocke et all 
Il99li iMaccacaro et all Il994|) . The EMSS appears to be 
euclidean and almost coincident with our curve for fluxes 
above 10 -13 erg cm -2 s _1 . 

We have also extended t he comparison to include 
samples at at fainter fluxes. iGeorgakakis et alJ l)2004a|) 
have computed the logN-logS of sources in the XMM - 
Newton/2dF survey, obtained with the EPIC instrument 
on board the XMM - Newton satellite. This survey covers 
an area of about 2.5 deg 2 to the flux limit ~ 10~ 14 erg 



m the 0.5 - 8 keV band (or F . 5 _ 2 ~ 5 x 10" 
erg cm" 2 s _1 ). They find two "normal" galaxies in their 
sample, that imply a density of < 1 source at Fx ~ 10~ 14 
erg cm" 2 s _1 , lower than what we find. 

However, the sample of lGeorgakakis et all 1 2004ah con- 
tains three additional galaxies (at z < 0.1 with Lo.5_s ~ 
10 42 erg s _1 ) that the authors do not consider because 
they might contain Low Luminosity AGNs (LLAGNs,see 
their Fig. 2). Since we cannot exclude that our sample 
also contains a few LLAGNs (see above and sectional), 
we should consider these objects for a better comparison 
with our sample. With the inclusion of these objects the 
logN-logS better matches ours. 

Also shown in Fig. are the constraints from the 
stacking analysis results of IGeorga kakis et all l|2003h at 
fainter fluxes, computed from optically selected galaxies 
at a mean redshift of z ~ 0.1. The point derived from the 
total sample considered is in excellent agreement with our 
logN-logS. 

Recently IGeorgakakis et all l|2004bh have presented a 
pilot sample of normal galaxies serendipitously detected 
in XMM - Newton public observations. They found 11 
"normal" galaxy candidates with luminosities below 10 42 
erg s _1 over an area of ~ 4.5 deg 2 . They find that the 
LogN-LogS derived from this sample (plotted in Fig. EJ) 
is again almost euclidean in slope, although at a smaller 
normalization than ours. 

We also plot the results from deep e r surveys, us 



ing da ta from lHornschemeier et all l|2003l). iNorman et al 
iported bv iRanalli et all l(2004h . iBauer et al 



as rer 



1 

20041) . derived from the Chandra Deep Fields. All these 
relations fall close to, although in general below, the ex- 
trapolation from our sample. Different authors derive dif- 
ferent slopes for their samples, but they are all consistent 
with the euclidean one. 

The lower normalizations found in these latter sam- 
ples could be explained in part by the combined effects 
of more stringent criteria to minimize contamination from 
the AGNs, even though of low luminosity, and of different 
relative occurence of the galaxy types (spiral/starburst vs 
early types). 

As already discussed, some residual contamination 
from low luminosity AGNs could be present in the sample 
we have considered, since we have little information on the 
optical spectra, a nd could only reject know n AGNs. We 
note however that lZimmerman n e^ aljl |200ll) have applie d 
the same criterium adopted bv lGeorgakakis et al.1 l|2004ah . 
namely an X-ray to optical luminosity ratio smaller than 
10 -2 , and that the EMSS sample, which is well studied 
optically, should not be contaminated by AGNs. We have 
nevertheless considered discarding sources in our sample 
that have a log(Fo.5-2/Fs) < —2. We have 7 objects that 
violate this limit, mostly at the high flux end. The re- 
sulting logN-logS relation is slightly steeper, but consis- 
tent with that presented in Fig. |5j and would not signif- 
icantly lower the normalization of the "euclidean" curve 
plotted. However the loca tion of the different points from 
IGeorgakakis et all 1 2004al) that consider/discard possible 
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contamination from LLAGNs give an idea of the possible 
uncertainties involved. 

The effect of different relative contributions from the 
early/late t ypes is more complicated to assess. The lowest 
flux points ijHornschemeier et all l2003j) are derived from 
late type galaxies, so th ey could underrepr esent the to- 
tal population. However. iBauer et all j20O4) suggest that 
the early type galaxies follow a flatter distribution, and 
in any case their LogN-logS is always below those from 
starburst/quiescent galaxies (see their Fig. 9), so their 
contribution could in fact be negligeable at the fainte r 
fluxes. In the range covered bv lGeoreakakis et all l|2004b|) . 
there is only one early type galaxy (but the sample is very 
small), while in our sample we have a sizeable fraction of 
early types (1/3 among the objects with a morphologi- 
cal classification), and the percentag e increases to ~50% 
in the "candidate normal galaxies" of lZimmermann et all 
l)200l[) . for which moreover there is no apparent signifi- 
cant different in the two slopes either. The stacking anal- 
ysis results from iGeorgakakis et alJ l|2008^ derived sepa- 
rately for E/SO and Sa-Scd bracket the extrapolation of 
the euclidean logN-logS obtained from brighter samples, 
with Sa-Scd in better agreement with the samples at low 
fluxes. Since the emission from early and l ate type galaxies 



is due to significantly d ifferent pro c esses_ 
Fabbiano et all Il992t iKim et all Il992 " 



Eskridge et al 



1995t IShaplev et all l2001f) , a different evolution is not out 
of the question. An assessment of the local LogN-logS sep- 
arately for the different morphological types would be a 
step forward to a better understanding of the properties 
of galaxies as a class, and would provide stronger con- 
straints for the investigation of normal galaxies at higher 
redshifts. Investigating this aspect is beyond the scope and 
the potential of this work. For the present time, we simply 
notice how remarkable it is that, in spite of the different 
selection criteria and instruments used to define all the 
samples considered, the surface densities of normal galax- 
ies is consistent with a single euclidean distribution for 
about 6 decades in flux (from ~ 10~ n to ~ 3 x 10~ 17 erg 
cm -2 s _1 ). 



8. Conclusions 

We present the results for an "almost serendipitous" sam- 
ple of 143 X-ray emitting normal galaxies selected from the 
cross-correlation of the BMW - HRI Catalogue and the 
LEDA database. Isointensity X-ray contours are overlayed 
onto the optical images for all galaxies and are presented 
in an atlas in Fig. [3 The X-ray characteristics of the sam- 
ple, listed in Tabled are derived in a uniform way and are 
used in comparison with other samples in the literature. 
We find that the general properties of the total sample are 
in good agreement with those already known for normal 
galaxies. 

We have also selected a complete subsample of 32 truly 
serendipitous sources in the local universe (z < 0.07), for 
which we derived the logN-logS distribution in the flux 
range between ~ 1.1 and 110 x 10~ 14 erg cm~ 2 s _1 , in 



the 0.5 — 2 keV energy band. We find that the relation is 
consistent with the euclidean distribution. 

Moreover, we find a good agreement between our logN- 
logS and those derived from ROSAT PSPC and Einstein 
data at similar or brighter fluxes and from XMM - Newton 
and Chandra at fainter fluxes: the overall distribution ap- 
pears to be consistent with an euclidean slope for about 
6 decades in flux, from ~ 3 x 10~ 17 to 10 -11 erg cm -2 
s _1 . The normalizations of different samples are consis- 
tent within a factor of ~ 2. 

Although with limited statistics, this work provides a 
first estimate of the number density of sources identified 
with normal galaxies in the nearby universe, in a flux range 
(both optical and X-ray) easily accessible for follow up 
detailed observations. This will allow us to provide a solid 
basis that will help in studying and classifying objects 
found in deeper surveys. 

While current efforts are mainly focused on probing 
the distant universe to determine the relevance of nor- 
mal galaxies as a class at very faint fluxes, the success of 
these studies depends also on the constraints given by the 
bright flux end of the number counts, which is still poorly 
studied. This is particularly relevant at the lowest fluxes 
where the number counts approach the total source den- 
sity (e.g. the detec t ions a nd fluctuation analysis results o f 
iMivaii k Griffiths! l)2002t . see lHornschemeier et all l|2002t 
and Fig. 7). The sample here derived will be therefore in- 
strumental to studies of both the galaxy cosmological evo- 
lution and the contribution of this class of sources to the 
X-ray background. 
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Fig. 2. First page of the figure available in electronic form, that presents the X-ray contours from smoothed X-ray 
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Gaussian function of a = 5" for point-like sources and of a — 10" for extended sources. 
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Appendix A: Results of optical observations 

The log of observations made at the Loiano Telescope in 
order to measure redshifts and magnitudes for some galax- 
ies in our sample (see section |SJ) is listed in Tabic HOI 

The seeing was about 2", so we used a long slit of width 
2.5" for spectra. The spectral range was 3800 A - 8700 A. 

For photometry we used the B - V - R filters of Johnson 
- Kron - Cousin system. 

We made use of the Image Reduction and Analysis 
Facility (IRAF) 6 for spectroscopic and photometric data 
reduction. 

We used all lines visible in the spectra to measure z, 
which is calculated from the average value. The redshifts 
obtained are reported in Table lA~2l Only absorption lines, 
typical of early- type galaxies, and no sign of nuclear ac- 
tivity were detected in the galaxy spectra. 

The atmosferic conditions were inadequate for photo- 
metric observations, however we list our raw estimate of 
magnitudes in Table I A. 31 and we use the value obtained 
for PGC 69291 in Tabled where it is reported in brackets. 

PGC 196693 
After obtaining the spectrum, we discovered that NED 
lists the galaxy SDSS J234817.99+010617.0 at the same 
position and with the same redshift we obtained, although 
it does not give a PGC name. We therefore identify 
PGC 196693 as SDSS J234817.99+010617.0. 

PGC 69291 

The galaxy is identified as NGC 7325 by LEDA and 
as NGC 7327 by NED, with a note suggesting that it 
could be a star. Our observations show 2 objects, a star 
and a galaxy at z — 0.0290 ± 0.0009 at the position of 
PGC 69291. Most of the optical flux is due to the star, as 
listed in Table lA~3l 

The uncertainty in the X-ray position (see also the 
map in Fig. FJJ does not allow us to identify the X-ray 
source with either object. Unfortunately the bad quality 
of our optical does not allow us to use the X- r ay - t o - 
optical ratio (see for example iMaccacaro et all Il988|) to 
discriminate between them, because it is consistent with 
both the star and the galaxy. We list the galaxy in the 
total sample but due to its faint estimated magnitude it 
is not included in the complete serendipitous sample. 

Appendix B: Table of total sample 

We report here the complete table with our total sample 
that will be available in the electronic form. For expla- 
nation of the columns refer to the sample Table ^ in the 
main text. 



6 IRAF is a general purpose software system for the re- 
duction and analysis of astronomical data. IRAF is writ- 
ten and supported by the IRAF programming group at the 
National Optical Astronomy Observatories (NOAO) in Tucson, 
Arizona. NOAO is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc. under cooperative 
agreement with the National Science Foundation. 
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Table A.l. Log of optical observations. 



Name 


Redshift 


PGC 65435 
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± 


0.0006 


PGC 69291 
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0.0009 
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Table A. 2. Average value of redshifts measured for the three galaxies observed. 
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Table A. 3. Estimated magnitudes for two of the three galaxies observed and for the star found near PGC 69291 
resulting from our observations with the Loiano Telescope. 
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